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ABSTRACT 

We present observations of the Lynds' dark nebula LDN 1111 made at microwave 
frequencies between 14.6 and 17.2 GHz with the Arcminute Microkelvin Imager (AMI). 
We find emission in this frequency band in excess of a thermal free-free spectrum 
extrapolated from data at 1.4 GHz with matched ui>-coverage. This excess is > 15a 
above the predicted emission. We fit the measured spectrum using the spinning dust 
model of Drain & Lazarian (1998a) and find the best fitting model parameters agree 
well with those derived from Scuba data for this object by Visser et al. (2001). 
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1 INTRODUCTION 



Recent pointed observations (Finkbeiner et al. 2002, 2004; 
Casassus et al. 2004, 2006; Watson et al. 2005; Scaife et al. 
2007; Dickinson et al. 2007) have provided some evidence 
for the anomalous microwave emission commonly ascribed 
to spinning dust (Drain & Lazarian 1998a,b). Although this 
emission was originally seen as a large scale phenomenon 
in CMB observations (see e.g. Kogut et al. 1996) it has 
been suggested that the emission occurs in a number of dis- 
tinct astronomical objects, such as dark clouds, and Hn and 
photo-dissociation regions. It often appears to be correlated 
with thermal dust emission as supported by the pointed ob- 
servations mentioned previously but it must be stated that 
this is not always the case, see for example Casassus et al. 
(2008). 

In spite of these predictions the evidence for anoma- 
lous microwave emission in compact objects is often contra- 
dictory. Early observations below 10 GHz of the molecular 
cloud LPH 96 (Finkbeiner et al. 2002), which showed a rising 
spectrum were later contradicted by observations at 31 and 
33 GHz which found emission consistent with an optically 
thin free-free spectrum extrapolated from lower frequencies 
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(Dickinson et al. 2006; Scaife et al. 2007). Although some 
evidence for an excess was found in a sample of Southern 
Hn regions (Dickinson et al. 2007) and more significantly in 
RCW175 (Dickinson et al. 2008), no emission inconsistent 
with free-free was found in a sample of Northern Hn re- 
gions (Scaife et al. 2008). Casassus et al. (2004) proposed a 
flux density of approximately 1 Jy from the Helix planetary 
nebula at 31 GHz to be in excess of a free-free spectrum ex- 
trapolated from lower frequencies. However, based on flux 
densities from the literature at 1.4, 2.7, 6.63 GHz (Ehman, 
Dixon, & Kraus 1970; Higgs 1971; Wall, Wright, & Bolton 
1976) and the recent WMAP 5-year densities at 23-94 GHz 
(Wright et al. 2008), which are all also ~ Uy, we suggest 
that there is no evidence for this reported excess in the flux 
density spectrum, although the nebula may be anomalous 
in other ways. 

In this letter we present observations of the Lynds' dark 
nebula LDN 1111, taken from the AMI sample of compact 
Galactic star formation regions (Scaife et al., in prep). This 
sample was selected from the SCUBA sample of compact 
Lynds' clouds (Visser, Richer & Chandler 2001). 

The spectra of dark clouds at gigahertz frequencies is 
poorly documented in all but a few cases (Finkbeiner et al 
2004; Casassus et al. 2006; Casassus et al. 2008). In those 
cases where cm-wave data is available (Casassus et al. 2006; 
Casassus et al. 2008) the behaviour of these objects has been 
found to be anomalous in a number of ways and in the case 
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of LDN 1622 (Casassus et al. 2006) to show a distinct excess 
of microwave emission. 

LDN 1111 (a = 21 h 40 m 30 s ,<5 = +57°48'00" J2000) lies 
on the inside edge of the heel of the large (w 3 degree) horse- 
shoe shaped Hn region IC 1396 (Sh 2-131; Sharpless 1959). It 
has no known IRAS association (Parker 1988) but has been 
studied in the submillimetre at 850 /im with the SCUBA 
instrument (Visser, Richer & Chandler 2001). An opacity 
class 6 object (A v > 5 mag; Lynds 1962), it is one of the 
most opaque dark nebulae. 



2 OBSERVATIONS 

2.1 Calibration and Data Reduction 

The AMI Small Array (SA) is a radio interferometer 
which observes in eight frequency channels in the band 
12-18 GHz at the Mullard Radio Astronomy Observatory, 
Lord's Bridge, Cambridge, UK. In practice, the lowest two 
frequency channels are generally unused due to a low re- 
sponse in this frequency range, and interference from geo- 
stationary satellites. AMI Consortium: Zwart et al. (2008) 
discusses the telescope in more detail. 

Observations of LDN 1111 were made in 13 hours over 
two days in 2007 November. Data reduction was performed 
using the local software tool reduce. This applies both au- 
tomatic and manual flags for interference and shadowing 
and hardware errors. It also applies phase and amplitude 
calibrations; it then Fourier transforms the correlator data 
to synthesize the frequency channels before output to disk 
in uv FITS format suitable for imaging in MPS. 

Flux calibration was performed using short observations 
of 3C286 near the beginning and end of each run. We as- 
sumed I+Q flux densities for this source in the AMI SA 
channels consistent with Baars et al. (1997) ~ 3.3 Jy at 
16 GHz. As Baars et al. (1977) measure I and AMI SA mea- 
sures I+Q, these flux densities include corrections for the 
polarisation of the source derived by interpolating from VLA 
5, 8 and 22 GHz observations. A correction is also made for 
the changing intervening air mass over the observation. From 
other measurements, we find the flux calibration is accurate 
to better than 5 per cent (Scaife et al. 2008; Hurley- Walker 
et al., in press). 

The phase was calibrated using hourly interleaved 
observations of the point source J2201+508 (a — 
22 h 01 m 43 s .5,<5 = 50°48'56".4), which has a flux density of 
0.3 Jy at 16 GHz. It was selected from the Jodrell Bank VLA 
Survey (JVAS; Patnaik et al. 1992). After calibration, the 
phase is generally stable to 5° for channels 4-7, and 10° for 
channels 3 and 8. In this work we use only channels 4-7 due 
to their superior phase stability. 

The FWHM of the primary beam of the AMI SA is 
» 20' at 16 GHz. The FWHM of the synthesised beam of the 
combined channel map towards LDN 1111 is 2.4' x 2.1' using 
natural weighting. 

2.2 Imaging 

The reduced visibilities were imaged using the AIPS data 
package. Dirty images were deconvolved using IMAGR which 
applies a differential primary beam correction to the CLEAN 
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Figure 1. Cleaned combined channel map of LDN 1111 at 
16 GHz. Contours arc increments of 1 mjy beam -1 , with the first 
contour at 1 mjy beam -1 . The primary beam FWHM of the AMI 
SAis shown as a black circle and the synthesized beam as a grey 
ellipse in the bottom left-hand corner. 
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Figure 2. Cleaned map of the uv sampled CGPS data at 1.4 GHz. 
Contours are increments of 1 mjy beam -1 , with the first contour 
at lmjy beam -1 . The primary beam FWHM of the AMI SA is 
shown as a black circle and the synthesized beam as a grey ellipse 
in the bottom left hand corner. 



components to account for the different frequency channels 
of the AMI instrument. Maps were made from both the com- 
bined channel set, shown here, and for individual channels. 

Two objects are labelled in both Fig. 1 and Fig. 2: ob- 
ject A, which is the dark cloud; and object B, which is the 
NVSS radio source J213955+574859. 
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Table 1. Flux densities for J213955+574859 



Frequency 


Flux density 


Reference 


(GHz) 


(mjy) 




0.327 


20.0±4.8 


WENSS; Rengelink et al. 1997 


0.408 


13.4±1.3 


this work; data from 






Taylor et al. 2003 


1.420 


12.2±1.25 


this work; data from 






Taylor et al. 2003 


1.420 


11.1±0.5 


NVSS; Condon et al. 1996 


14.631 


3.48±0.56 


this work 


15.381 


2.88±0.60 


this work 


16.130 


3.02±0.55 


this work 


17.150 


2.93±0.57 


this work 



2.3 Radio Spectra 

Since the frequency coverage of AMI tells us nothing about 
the spectral behaviour at longer radio wavelengths we must 
combine our new data with existing archival data. Ideally 
total power measurements are required, with the same or 
better angular resolution than AMI. Data of this type al- 
low sampling in the uv plane to match exactly the mea- 
sured angular scales of AMI. Here we use data from the 
Canadian Galactic Plane Survey (hereinafter CGPS; Taylor 
et al. 2003) at 1.42 GHz. These data are a combination of 
single-dish and synthesis data, and provide a total power 
measurement of this region with a resolution of « 1'. In 
order to mimic an AMI observation, these data are modu- 
lated by the primary beam of AMI before sampling in the 
uv plane to match the AMI uv coverage. These data arc 
then processed in the same way as the true AMI visibilities 
to give the final sampled image. Fig. 2 shows the sampled 
data which result from matching the visibility coverage to 
AMI channel 4. 

Spectra of both LDN 1111 and the radio source 
J213955+574859 are presented, identified as objects A and 
B respectively in Fig. 1. Errors on the data points in these 
spectra are calculated using a contribution for the rms noise 
in each channel, a conservative 5% flux calibration error and 
a contribution for the uncertainty in flux extraction. The 
morphology of LDN 1111 and its surroundings makes this ex- 
traction subject to the fitting area used. To account for this 
in the errors 10 independent 'fits' are made. The recorded 
flux density is the mean of these measurements and the 
variance is referred to as crf t . These errors are combined in 
quadrature: as = \/cr?ms + (0.05Si) 2 + o\ t . The flux extrac- 
tion is performed using the fitflux program (Green 2007), 
which takes a user defined polygonal aperture and performs 
aperture photometry after subtracting a twisted plane back- 
ground level. These errors are dominated by the 5% calibra- 
tion uncertainty, with the map r.m.s. being ~0.2mJy per 
channel. 

The spectrum of J213955+574859 is presented in or- 
der to demonstrate the flux calibration of AMI, see Fig. 3. 
The data points for this plot are given in Table 1. In ad- 
dition to data points from the literature, flux densities at 
408 MHz and 1.42 GHz from archival CGPS data are also 
plotted. These data have been extracted using the fitflux 
program as described above. Spectral indices are calculated 
using firstly the AMI data points plus the literature data, 



q = 0.52; and secondly using all the data points, a = 0.46. 
These fits are shown as dot-dash and dashed lines respec- 
tively. Both are consistent with a^Ss ~ 0-43 ± 0.18. 

For J213955+574859, which is approximately point- 
like, we expect a negligible amount of flux loss from ex- 
tended structure across the AMI bandwidth. In the case 
of LDN 1111 the same is not true. Using the CGPS map 
at 1.42 GHz we find > 50% flux loss: the integrated flux 
from the original unsampled map, calculated as above, is 
Si. 42, total ~ 25.4 ± 2.7 mjy. However, since the morphology 
of the source varies between 1.4 and 16 GHz we have not 
used the CGPS data as a model to calculate flux loss. 

Flux densities for LDN 1111 are listed in Table 2. It is 
immediately obvious that the flux densities for LDN 1111 
between 14 and 18 GHz are inconsistent with an optically 
thin free-free spectrum extrapolated from 1.42 GHz. At the 
high frequency end of the AMI spectrum there is an excess 
of almost 50mJy, > 15<rs (« 250u rms ). 

Although an absolute measure of the flux loss is not 
feasible, it is possible to make an estimate of how the rel- 
ative flux loss will affect the shape of the microwave spec- 
trum. This is done using a multi-variate Gaussian model 
for the source with dimensions of the deconvolved source at 
14.63 GHz (AMI channel 4) as given by the aips task jm- 
fit, 7' x 4'. This model is then sampled in uv using the 
exact coverage of the individual channels, and their flux loss 
relative to the lowest channel is recorded. We assume an 
uncertainty of 10% on the percentage flux loss and propa- 
gate the errors. The results of this calculation are shown in 
Fig. 4 and tabulated in Table 2. It should be noted that this 
is a conservative uncertainty and the relative flux loss will 
only be in error should the morphology change significantly 
between 14 and 17 GHz. 

Fig. 4 includes a data point at 353 GHz (850 pm) from 
the Scuba instrument (Visser, Richer & Chandler 2001) for 
illustartive purposes. This data point is uncorrected for flux 
loss. These Scuba data at 850 //m are chopped at 2.5 arcmin 
and consequently we would assume that information on 
scales larger than this is lost. To correctly estimate the flux 
loss at this frequency we would require information on the 
morphology of the vibrational dust emission. Parker (1988), 
from whose paper the Scuba sample were selected, lists the 
dimensions of LDN 1111 as 1.7 x 1.1 arcmin. A multivariate 
Gaussian model based on this data would suggest a flux loss 
of approximately 73%. 

Source extraction methods applied to the WMAP data 
at 23-91 GHz (Lopez-Caniego et al. 2007) are only able 
to give us upper limits on any compact emission towards 
LDN 1111, with the lowest limit being at 41 GHz, Sai < 
0.30 Jy. This value is consistent with the noise level at that 
frequency. 



3 ANALYSIS 

The emission seen between 14.63 and 17.15 GHz by AMI 
is clearly in excess of a simple free-free spectrum extrapo- 
lated from the CGPS data at 1.4 GHz. A number of possi- 
bilities may throw some light on the nature of this excess. 
The thermal (vibrational) dust spectrum of proto-planetary 
disks around T-Tauri stars is expected to extend into the cm- 
regime. In these circumstances grain growth has increased 
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Figure 3. Spectrum of the source J213955+574859. Data points 
from the literature arc shown as crosses and are described in the 
text. Data points extracted from the CGPS data are shown as 
un-fillcd triangles. Data points from the AMI are shown as filled 
circles. 



freq/GHz 

Figure 4. LDN1111: Data points are listed in Table 2, with 
the addition of the SCUBA flux density at 850 (im. A free-free 
spectrum as described in the text and a grcybody spectrum scaled 
to the Scuba data point are shown as dashed lines. The best-fit 
spinning dust model of Draine & Lazarian (1998) is shown as a 
dot-dash line and the total spectrum as a solid line. 



Table 2. Flux densities for LDN1111 



Frequency Flux density Relative Corrected 
flux loss flux density 
(GHz) (mjy) (%) (mjy) 



1.420° 


7.2 ±0.7 




7.2±0.7 


14.631 


50.7±2.9 




50.7±2.9 


15.381 


53.9±3.2 


10 


59.9 ±6.8 


16.130 


54.5±3.3 


17 


65.7 ±7.6 


17.150 


55.5±3.3 


24 


73.0 ±8.4 



a uv sampled CGPS data point 



the population of cm-sized grains, or pebbles, causing (5 to 
approach zero and the greybody spectrum to fall off with 
an index approaching 2. However the size of these disks is 
small and would be unlikely to be resolved by either AMI or 
Scuba. LDN1111 is quite obviously elongated in the AMI 
data, as it is in the Scuba map which has a resolution of 
14arcsec. Although we do not deny the possibility of there 
being such a disk embedded within LDN 1111 or projected 
along the line of sight through this object, in the absence of 
further evidence we choose to pursue an alternative expla- 
nation. 

A further possibility for a superposition of LDN 1111 
and another object is that there may be an Hn region, with 
a density such that it is still optically thick at 16 GHz, whose 
flux is contributing to the spectrum. An Hn region with a 
turn-over frequency > 30 GHz, as would be required here, 
would possess an emission measure in excess of 4 x 10 pc 
cm~ 6 . This would put it in the regime of ultra-compact and 
hyper-compact Hn regions. The spectrum of these objects 
is inverted with an index of ~ — 2. The AMI data points 
alone have a spectral index of a = —2.86 ± 0.35, rather 
steeper than would be expected. Indeed hyper-compact Hn 
regions exhibit slightly shallower spectra (Franco et al. 2000) 
which extends even to mm-wavelengths. The lack of a visible 
source in the WMAP data towards LDN 1111 would suggest 
that this is not the case here. In addition hyper-compact 
Hn regions tend to have very broad radio recombination 
line (RRL) widths, typically in excess of 50 km s _1 (Kurtz 
2005). Ultra-compact Hn regions also show broadened line 



widths, of typically 30-40 km s - . The RRL measurement of 
Heiles, Reach & Koo (1996) which encompasses LDN 1111 
has a width of only 21.4 km s _1 , although we note that 
this measurement is not directed at LDN 1111 specifically. 
An emission measure of this magnitude would also imply 
an Hn region mass of around 400 M . The total mass of 
the cloud from the Scuba data (Visser et al. 2002) is 0.3 M 
only, which given the calculated flux loss in this data implies 
an upper limit on the mass of 1 M . 

A third possibility for the excess emission seen in the 
AMI data is that of dipole emission from rapidly rotating 
small dust grains (Drain & Lazarian 1998a,b). The obser- 
vational evidence for this emission mechanism has been ex- 
plored in the introduction to this letter. It is a relatively 
new emission mechanism and the evidence for its existence 
is not conclusive. We assess the possibility that the emission 
seen in LDN 1111 arises as a consequence of spinning dust 
by comparing the data to the model of Draine & Lazarian 
(1998a;b), which has been used extensively in the past. To 
make this comparison we use the MCMC based software 
metro (Hobson & Baldwin 2004) to find the best-fit pa- 
rameters. 

We fitted a model which has a free-free component nor- 
malized to give the flux density at 1.4 GHz from the sampled 
CGPS data with a spectral index of a = 0.1. To this we add 
a spinning dust component scaled from the DL98 molecular 
cloud model. 

The model is parameterized by the column density, 
iV(H2), and the angular size of the object. Visser et al. 
(2001) calculate the average and peak column densities for 
LDN 1111 as 5 and 13xl0 21 cm~ 2 , respectively. The results 
of our model fitting are N(H 2 ) = 6.72±0.58xl0 21 cm~ 2 , and 
A8 = 5.38±0.26 arcmin. These values provide a \ 2 of 1.03 
(79% probability). However, given the small number of de- 
grees of freedom in this case we do not propose this statistic 
as being conclusive. The resulting model is shown in Fig. 4. 
The column densities agree well with those of Visser et al., 
although as might be expected there is a degeneracy between 
the two parameters. 
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4 DISCUSSION AND CONCLUSIONS 

Given the available evidence we propose the excess of emis- 
sion we see towards the dark cloud LDN 1111 in the mi- 
crowave band to be a result of emission from small spin- 
ning dust grains. This excess is quantified relative to a flux 
density derived from lower frequency data sampled to give 
the same uv coverage and consequently measuring the same 
angular scales on the sky. Lacking further data at lower fre- 
quencies with suitable angular resolution this one point is 
the basis for the assumption that there is an excess in the 
AMI frequency band. We therefore qualify this excess by 
using the original CGPS data to measure the total power, 
i.e. power on all scales, flux density at 1.42 GHz to be 
Si. 42, total = 25.4 ± 2.7 mjy. Since the object is extended 
this value may be regarded as an upper limit on the mea- 
sureable flux at 1.42 GHz. In addition the total power CGPS 
data at 408 MHz may then be used to constrain better the 
lower frequency spectrum: £408, total = 21.6 ± 11.8 mjy. The 
large error on this flux density is dominated by the fitting 
uncertainty. In combination these two flux densities yield a 
spectral index of a = —0.16 ± 0.12, which predicts a total 
power flux at 16 GHz of Sw, total = 36.9tlo'.i mjy. This sug- 
gests that even in terms of total power there is a clear excess 
at 16 GHz. 

Relative to a canonical free-free spectrum extrapolated 
from our sampled CGPS data at 1.42 GHz we see an excess 
of Sexccss = 60.1 ± 7.6 mjy, « 8a, assuming a fixed spectral 
index of a = 0.1. Using the spectral index derived from 
the CGPS total power maps at 408 MHz and 1.42 GHz then 
there is an excess of S e ^ cess r = 55.10 mjy, which is still > 7a. 

In conclusion, we have presented observations of the 
Lynds' dark nebula LDN 1111 at frequencies of 14.6- 
17.2 GHz. These measurements show an excess of emission 
towards this object relative to an extrapolated free-free 
spectrum at a significance of w 15<7. We have proposed that 
this excess may be due to emission from small spinning dust 
grains and find that it is well-described by the model of 
Draine & Lazarian (1998a;b). 
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